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SECTION  I 


INTRODUCTION 

Accurate  knowledge  of  accelerometer  nonlincarities  is  required  in  the 
application  of  guidance  control  to  today's  missiles  and  space  vehicles.  Though  some 
attempts  have  been  made  to  determine  these  coefficients  on  a  precision  linear  shake 
table,  the  only  successful  method,  thus  far,  has  been  by  use  of  a  precision  centrifuge. 
However,  the  use  of  a  precision  centrifuge  does  not  guarantee  that  the  nonlinear 
coefficients  will  be  accurately  determined  since  there  are  several  significant  error 
sources  that  must  be  controlled. 

A  data  reduction  procedure  to  minimize  the  effects  of  misalignments  and  radius 
arm  uncertainties  has  been  described  by  Evans  and  Fuhrman  (Ref.  1)  and  is  widely 
used  by  the  industry.  In  this  procedure,  quadratic  equations  are  separately  fitted 
to  accelerometer  output  data  taken  in  two  positions  (acceleration  vector  along  plus 
and  minus  input  axis)  by  the  method  of  least  squares.  The  constant  and  linear  terms, 
so  obtained,  are  subtracted  from  the  data  and  a  best  fit  cubic  equation  is  fitted  to 
the  resulting  combined  residuals, 

Though  the  above  procedure  is  mathematically  incorrect,  it  works  fairly  well 
for  determining  the  nonlinear  coefficients  because  these  coefficients  are  relatively 
small  and  do  not  need  to  be  known  with  great  precision. 

Evans  and  Fuhrman  have  also  outlined  an  iteration  procedure  to  improve  the 
accuracy  with  which  the  nonlinear  coefficients  may  be  obtained. 

In  this  paper,  a  simpler,  more  direct,  more  accurate,  and  mathematically 
correct  data  reduction-procedure  is  presented.  It  is  a  procedure  in  which  all  the 
data  (bipolar)  are  combined  at  cnce  for  the  optimum  separation  of  even  and  odd  terms. 
It  does  involve  a  larger  matrix  of  simultaneous  equations  than  the  method  of  Evans 
and  Fuhrman  but  with  modem  computers,  even  desk  computers,  this  is  no  great 
problem. 
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The  data  reduction  procedure  has  been  applied  not  only  for  obtaining  the  input 
axis  nonlinear  coefficients  but  also  for  obtaining  all  th<>  other  nonlinear  coefficients 
of  the  assumed  model  equation.  Kq.  (2-1).  Methods  for  establishing  the  uncertainties 
of  the  coefficients  and  the  adequacy  of  the  model  equation  for  particular  applications 
are  also  presented. 

U  is  shown  that  the  radial  acceleration  gradient,  inherent  to  a  centrifuge,  has 
no  effect  on  a  rigid  body  pendulum  but  that  the  rotation  may  affect  the  results  due  to 
product  of  inertia  torques.  In  addition,  there  is  a  discussion  on  the  severe  limitations 
imposed  on  the  determination  of  these  coefficients  by: 

(1)  misalignment  angles 

(2)  uncertainties  in  the  length  of  the  radius  arm 

(3)  relative  magnitudes  of  the  coefficients 

(4)  angular  compliance  of  the  centrifuge  arm  and  of  the  mounting 
fixture  about  all  three  principal  axes. 

The  above  limitations  are  reasonably  independent  of  the  data  reduction  procedures. 
Limitations  (1)  and  (2),  above,  were  discussed  by  Evans  and  Fuhrman  wtth  regard  to 
obtaining  bias  and  scale  factor  but  not  with  regard  to  their  effects  on  the  determination 
of  the  nonlinear  coefficients,  in  general. 

The  misalignments  include  not  only  the  misalignments  of  the  true  input  axis 
with  respect  to  the  accelerometer  mounting  surface(s)  but  also  the  misalignments  of 
fixture  mounting  surface (s)  and  the  centrifuge  mounting  surface(s)  with  respect  to  the 
centrifuge  and  gravity  acceleration  vectors.  In  addition,  misalignments  can  be  introduced 
by  tightening  of  mounting  screws,  by  dirt  on  the  mounting  surfaces,  or  by  ether  causes. 

Depending  on  the  magnitude  of  the  nonlinear  co-fficients,  it  was  found  that  the 
angular  compliance  of  the  centrifuge  arm  and  mounting  fixture  should  not  exceed  about 
1  or  2  are  second  per  g;  about  any  axis. 
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SECTION  II 


MODEL  EQUATION 


For  convenience  in  this  presentation,  the  model  equation  for  a  linear,  pendulou 
nongyroscopic,  torque  balanced  accelerometer  is  expressed  in  terms  of  the  applied 
acceleration  components  along  the  true  input,  pendulous,  and  output  axes  rather  than 
along  the  reference  axes  as  given  in  Ref.  2.  The  misalignments  6p  and  6q  (Ref.  2) 
of  the  true  input  axis  with  respect  to  the  input  reference  axis  are  included  later  in  the 
misalignments  of  the  accelerometer  axes  with  respect  to  the  centrifuge  and  gravity 
acceleration  vectors. 


K0  +  ai  *  K21af  +  K3iaf  + 


+  Ko„a?  +  K,  a®  +  K.a.a 
2o  o  3o  o  ip  1  p 


(2-1) 


where 


A 

E 


aFap’ao 


=  acceleration  indicated  by  the  accelerometer  -  g1 

-  accelerometer  output-output  units 

o 

=  applied  acceleration  components  along  the  true  input, 

pendulous,  and  output  (pivot)  or  flexure  axes,  respectively  —  g. 

-  bias— g 

=  scale  factor  —  output  units/g 

2 

=  second  order  input  axis  coefficient —g/g 


*g  is  a  unit  of  acceleration.  For  convenience,  it  is  usually  taken  to  be  the  local  value 
of  gravity  though  some  standard  value  may  be  chosen,  if  desired. 

“Applied  acceleration  refers  to  nongravitational  acceleration  since  an  accelerometer 
cannot  sense  the  acceleration  of  free  fall.  The  attractive  force  of  gravity  acting  on 
an  earthbound  accelerometer  is  equivalent  in  its  effect  to  an  upward  applied 
acceleration  of  one  g. 


K3i 

V  K2o 

K3p’  So 


K 


ip’  Kpo’  Koi 


K. 


O 

third  order  input  axis  coefficient  -g/g 

2 

second  order  cross-axis  coefficients  -  g/g 

third  order  cross-axis  coefficients  -  g/g 

2 

cross-coupling  coefficients  -  g/g 

2 

torquing  power  (heating)  coefficient  -  g/g 


The  output  units  may  be  in  volts,  mA,  pulses/sec  or  other  convenient  units. 

Linear  cross-axis  sensitivity  terms,  K  a  and  K  a  have  not  been  included  in  the 

p  p  oo 

model  equation  since  they  cannot  be  distinguished  from  misalignments  of  the  input 
axis  about  the  pendulous  and  the  output  axes. 


The  input  reference  axis  (IRA),  the  pendulous  reference  axis  (PRA),  ar.d  the 
output  reference  axis  (ORA)  form  a  cartesian,  orthogonal  coordinate  system  which  is 
generally  indicated  by  case  markings  and/or  mounting  surface(s).  The  true  input 
axis  (IA),  pendulous  axis  (PA),  and  output  axis  (OA),  to  which  Eq.  (2-1)  refers,  are 
generally  slightly  misaligned  with  respect  to  the  reference  axes  and  can  only  be 
determined  by  test.  The  output  axis  (OA)  is  mechanically  determined  by  the  flex  or 
pivot  axis  and  the  positive  direction  is  chosen  arbitrarily.  The  positive  direction  for 
PA  is  orthogonal  to  OA  and  in  the  direction  from  the  flex  or  pivot  axis  through  the 
centroid  of  the  proof  mass  when  the  proof  mass  is  at  its  null  position.  The  three 
axes  form  a  right-hand  orth ogonal,  coordinate  system  with  the  origin  at  the  centroid 
of  the  proof  mass  and  such  that 

OA  x  IA  =  PA  (2-2) 

where  OA,  IA,  and  PA  are  unit  vectors  along  OA,  IA,  and  PA,  respectively. 

If  the  torquing  power  is  kept  constant,  as  it  is  in  some  pulse  torqued  systems, 
then  the  coefficient  Kt  of  Eq.  (2-1)  may  be  zero.  Terms  may  be  added  to  or  deleted 
from  Eq.  (2-1)  as  appropriate  for  the  type  of  accelerometer  and  the  requirements  of 
the  application.  The  data  reduction  principles  outlined  herein  may  be  applied  to  any 
reasonable  model  equation. 


3. 


Fherc  may  be  other  sources  than  torquing  power  for  this  term. 
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SECTION  III 


ACCELERATION  GRADIENT  EFFECTS 

In  centrifuge  testing  of  an  accelerometer,  it  is  evident  that  the  acceleration 
field  is  nonuniform;  it  varies  linearly  with  the  radius  from  the  axis  of  rotation.  The 
question  is  frequently  asked:  What  effects,  if  arty,  does  this  linear  acceleration  gradient 
have  on  the  accelerometer? 

Construct  a  right-hand,  orthogonal  coordinate  system  XYZ  that  rotates  with 
the  centrifuge  as  shown  in  Fig.  1.  The  Z  axis  coincides  with  the  rotation  axis  and 
the  X  axis  passes  through  the  center  of  gravity,  C,  of  the  proof  mass  (the  proof  mass 
includes  all  pendulous  parts  and,  in  general,  must  be  treated  as  a  compound  pendulum). 


2 


Figure  1.  Proof  Mass 
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Construct  an  xyz  coordinate  system  with  the  center  of  gravity,  C,  of  the  proof 
mass  as  the  origin  and  parallel  to  the  XYZ  system  as  shown  in  Fig.  1. 

Let 


i,  j.k 


R 


M 


=  unit  vectors  along  X,  Y,  Z  axes,  respectively.  Also 
along  x,  y,  z  axes,  respectively 

=  radius  vector  from  the  Z  axis  to  the  center  of  gravity,  C, 
of  the  proof  mass. 

=  projection  on  XY  plane  of  the  radius  vector  from  the  Z  axis 
to  any  point  P  of  the  proof  mass 

=  radius  vector  from  point  C  to  point  P 

=  components  of  r  parallel  to  the  x,  y,  and  z  axes  respectively 

=  steady-state  angular  velocity  of  the  centrifuge  about  the 
Z  axis  with  respect  to  inertial  space 

=  mass  of  the  proof  mass. 


Consider  a  differential  mass,  dm,  at  point  P.  The  centripetal  force  acting  on 
the  differential  mass  is 

dF  =  -R2  R  dm  (3-1) 

But,  from  Fig.,  1,  it  is  evident  that 

R  =  RcI+  r  -rzk  (3-2) 

Therefore,  if  we  substitute  Eq.  (3-2)  in  Eq.  (3-1)  and  integrate  over  the  entire 
body  of  the  proof  mass,  we  find 


F  =  -fl2[lRc  /dm  +/  r  dm  -  k/r^  dm  J 

But / dm  =  M,  the  total  mass  of  the  proof  mass,  and  the  other  two  integrals  define 
the  center  of  gravity  of  the  proof  mass  and  are  equal  to  zero.  Therefore 


F  =  -a  R„ M  i 


(3-3) 
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and  the  acceleration  gradient  has  no  effect  on  the  central  force  vector  acting  at  the 
centroid  of  the  proof  mass. 

Let  us  now  determine  if  the  acceleration  gradient  causes  any  moments  to  act  on 
the  proof  mass.  The  moment  of  the  centripetai  force  acting  on  the  differential  mass 
about  the  centroid,  C,  is 

dT  =  rxdF  =  -fi2  dm  r  x  R  (3-4) 


Substitute  Eq.  (3-2)  in  Eq.  (3-4)  and  integrate  over  the  entire  body  of  the  proof  mass. 
t  =  -a2  [rJ?  x  i  dm  ♦A  xrdm«/rxr2  k  dm  j 


or 


T  =  -fl' 


He/' 
Tfry r*  dm  J 


r  dm  +  R  i  /  r  dm  +  T I  r  r  dm 
y  c  V  z  *  J  z  x 


(3-5) 


The  first  two  integrals  in  Eq.  (3-5)  define  the  centroid  and,  therefore,  are  equal 
to  zero.  The  next  two  integrals  define  products  of  inertia  of  the  proof  mass.  Therefore 

T  =  <V=r-'zxT»n2  I3-6) 

where  I  and  I  are  products  of  inertia  of  the  proof  mass.  This  torque  is  a 
yz 

dynamic  unbalance  torque  produced  by  rotation  about  the  Z  axis  or  an  axis  parallel 

to  the  Z  axis.  If  the  z  axis  is  a  principal  axis  of  the  proof  mass,  then  I  =  I  =0. 

yz  zx 

The  construction  of  most  accelerometers  is  such  that  IA,  PA,  and  OA  are 
nominally  principal  axes  of  Inertia  of  the  proof  mass.  In  all  mounting  positions 
considered  in  this  paper,  the  centrifuge  rotation  axis  is  nominally  parallel  to  IA, 

PA,  or  OA  and  no  product  of  inertia  effects  are  included.  In  general,  only  product 
of  inertia  torque  components  about  OA  need  be  considered  and  these  vary  linearly 
with  acceleration. 
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SECTION  IV 


ACCELERATION  COMPONENTS 


In  this  section,  the  effects  of  earth's  rate,  nonverticality  of  the  centrifuge  axis, 
angular  deflections  of  the  centrifuge  radius  arm,  and  misalignments  of  t.he  true 
accelerometer  axes  with  respect  to  the  accelerometer  reference  axes  are  considered. 
It  is  assumed  that  the  accelerometer  reference  axes  are  defined  with  respect  to  the 
centrifuge  and  gravity  acceleration  vectors. 

The  following  notation,  in  addition  to  that  previously  defined,  will  be  used: 


a  ,  a  ,  a 
x’  y*  z 


a  ,  a  ,  a 
u  v  w 


T 


=  centripetal  acceleration  at  centroid  of  proof  mass  —  g 

=  components  of  applied  acceleration  along  the  x,  y,  and 
z  axes,  respectively  —  g 

=  components  of  applied  acceleration  along  the  u,  v,  and 
waxes,  respectively  —  g 

=  components  of  applied  acceleration  along  IA,  PA,  and 
OA,  respectively 

o 

=  magnitude  of  g  -  ft/sec 

=  compliance  coefficients  of  mounting  fixture  and  centrifuge 
about  the  x,  y,  and  z  axes,  respectively  —  rad/g 

=  period  of  revolution  with  respect  to  earth-sec 


^x’  Yy’  Yz 


€ 


zv 
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=  angular  deflections  of  mounting  fixture  and  centrifuge 
about  the  x,  y,  and  z  axes,  respectively,  as  a  function 
of  ac  -  rad 

=  angle  between  Z  axis  and  local  vertical  —  rad 

=  angular  displacement  of  the  centrifuge  arm  with  respect 
to  an  earth-fixed  reference  plane  containing  the  Z  axis  —  rad 


Preceding  page  biaok 
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A 

9, 


u 


small  successive  Euler  angle  displacements  about  IRA, 
PRA,  and  ORA  to  transform  acceleration  components 
along  accelerometer  reference  axes  to  components  along 
the  true  axes  -  rad 

-  astronomic  latitude  at  centrifuge  -  deg 
=  angular  velocity  of  the  centrifuge  relative  to  earth-rad/sec 
=■  angular  velocity  of  earth  relative  to  inertial  space-rad/sec 


6 

The  earth  sidereal  rate  is  quite  small  (wc  -  72.9211  x  10  rad/sec)  and  may 
be  neglected  in  the  calibration  of  some  accelerometers,  except  on  centrifuges  with 
long  radius  arms  and  at  low  centrifuge  rates  of  rotation.  The  angular  velocity  cf 
the  centrifuge  with  respect  to  inertial  space  is 


9  =  9r  +  uc  sin  A 

and  the  centripetal  acceleration  at  the  proof  mass  centroid  is 

Rc^2  Rc(Rr  +we  sinX)2  „  RcRr2  /  2u>esinA^ 
ac  =  “  =  G  5  "  G  X  +  ~r  ) 


(4-1) 


Thus  the  vertical  component  of  earth  rate  at  the  middle  latitudes  would  result  in  an 
error  of  approximately  10-2/,/  a  %  for  R  =  32  ft  and  varies  approximately  as  the 

v*  v 

square  root  of  the  radius. 


The  centripetal  acceleration  is  a  function  of  the  angular  velocity  and  the  radius 
arm.  The  angular  velocity,  corrected  for  earth  rate,  if  necessary,  may  be  determined 
quite  precisely  by  measuring  its  period  of  rotation,  T,  over  several  revolutions.  The 
radius  to  the  centroid  of  the  proof  mass  is  more  difficult  to  obtain.  Some  centrifuges, 
such  as  the  100-inch  centrifuges  at  Holloman  AFB,  have  a  cage  rotating  about  a  vertical 
axis  which  permits  a  fairly  accurate  determination  of  the  radius  arm.  An  error  of 
0. 1  percent  in  the  radius  arm  will  result  in  an  error  of  approximately  0.2  percent  in 
a  quadratic  term  and  about  0,3  percent  in  a  cubic  term,  which  is  usually  quite 
acceptable.  The  centripetal  acceleration  at  the  proof  mass  centroid  is 
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(4-2) 


a 


c 


R 


+ 


Tea  sin  \ 

e _ 

v 


\ 

/ 


Due  to  rotational  stresses  that  vary  as  the  square  of  the  centrifuge  rate,  the 

radius  arm  length  may  change.  It  would  probably  be  safe  to  assume  that  arm  stretch, 

o 

whether  positive  or  negative,  is  proportional  to  Q  .  It  is  comparatively  easy  to 
accurately  measure  arm  stretch,  so  it  will  be  assumed  that  this  is  done  and  used  in 
the  data  reduction  procedure,  if  required.  In  general,  measurements  of  arm  stretch 
dc  not  include  changes  that  are  due  to  bending  of  the  mounting  fixture. 

The  centrifuge  axis  in  a  well-installed  and  adjusted  centrifuge  should  be  very 
close  to  vertical  so  that  e  «  1.  Assuming  that  this  angle  is  independent  of  n  , 
which  is  not  necessarily  true,  e  may  be  determined  by  observing  a  low  threshold, 

Z  V 

high  resolution,  precision  hubble  level  mounted  parallel  to  the  X  axis.  The  bubble 
level  reading  should  vary  sinusoidally  as  the  centrifuge  is  rotated  from  one  position 
to  another,  i.e. 


°b  =  °o  +  €zvCOs  ( 6  +  /3) 

where 

=  angular  reading  of  bubble  -  rad 

bo  =  offset  angle  (bias)  of  bubble  -  rad 

/3  =  phase  angle  with  respect  to  reference  plane  -  rad 

The  magnitude  of  is  equal  to  one-half  the  variation  in  over  a  complete 
revolution. 


As  mentioned  previously,  the  centrifuge  axis  of  rotation  may  vary  with  speed 
but  probably  not  by  significant  amounts.  The. acceleration  components  in  the  xyz 
system  are 


ax  “  "ac  •  SvCO8<0+0> 
a  «zv  sln  (»  +  P) 


a„  =  cos  e, 


zv 


(4-3) 
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From  Eqs.  (4-3),  it  is  seen  that  the  x  and  y  components  of  acceleration 
arising  from  nonverticality  of  the  centrifuge  axis  average  out  to  zero  over  an  integral 
number  of  revolutions.  However,  it  is  best  to  have  the  axis  vertical  within  one  or  two 
arc  minutes  so  that  the  output  reading  is  not  necessarily  tied  down  to  averaging  over  an 
integral  number  of  revolutions  and  any  rectification  effects  are  negligible.  It  will  now 
be  assumed  that  the  centrifuge  axis  is  adjusted  to  within  one  or  two  arc  minutes  of  the 
vertical.  The  average  acceleration  components  along  the  x,  y,  and  z  axes,  in  £  units, 
are: 


Due  to  bending  of  the  centrifuge  arm  and  of  the  mounting  fixture  induced  by 
dynamic  unbalance,  aerodynamic  effects,  or  heating  effects,  the  accelerometer 
mount  may  deflect  angularly  as  a  function  of  acceleration  level.  As  a  first  approxi¬ 
mation,  it  will  be  assumed  that  angular  deflections  will  be  proportional  to  the 
acceleration  level,  i.e. 


Tx  =  kx  ac 


yy  1  y  ac 
yz  ~  kz  ac 


(4-5) 


In  a  well  constructed,  dynamically  balanced  centrifuge  with  a  rigid  mounting 
fixture,  the  angular  compliance  should  be  less  than  two  arc  seconds  per  g.  For 
example,  the  compliance  of  the  MIT/CSDL  32-foot  arm  centrifuge  was  found  to  be 
less  than  0.2  arc  seconds  per  g.  The  compliance  of  the  Holloman  100-inch  arm 
centrifuge  about  the  y  (tangential)  axis  has  been  found  to  be  less  than  0.25  arc  seconds 
per  g  on  some  centrifuge  runs  but  as  much  as  4  arc  seconds  per  g  on  other  runs.  The 
latter  variations  may  be  due  to  variations  in  the  dynamic  unbalance  of  the  birdcage 
for  the  different  setups.  It  should  be  noted  that  the  contribution  of  the  mounting 
fixture  to  the  angular  compliance  may  vary  as  a  function  of  the  accelerometer 
mounting  position. 
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Construct  a  right-hand,  orthogonal  coordinate  system  uw  derived  from  the 

xyz  system  (see  Fig,  2)  by  the  Euler  angle  rotations  y  ,  y y  .  The  uv\v  system 

x  j  z 

coincides  with  the  xyz  system  when  a  =  0. 


Figure  2.  Coordinate  Systems 

Though  angular  displacements  about  orthogonal  axes  are  not  commutative,  they  may 
be  treated  as  commutative  when  they  are  very  small,  say  less  than  five  arc  minutes. 
Using  small  angle  approximations,  the  acceleration  components  in  the  uvw  system 
are: 
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Next  use  small  successive  Euler  angle  displacements  ©.,  ©p,  ©0  to  transform 
acceleration  components  along  the  accelerometer  reference  axes  to  acceleration 
components  along  the  true  accelerometer  axes.  The  resulting  relationships  are 
given  in  Eqs.  (4-8). 

a. -A.  coso  coso  +A  (sin  o.  sin©  cos  ©„  +  sin  ©„  cos  ©.) 

11  p  op  i  p  o  o  1 

A  (sin  o  sin  o.  -  cos  o.  sin  o  cos  ©1 

o  O  1  l  p  o 

a  Ap  (cos  0Q  cos  Oj  -  sin  o.  sin  op  sin  oq)  +  Aq  (cos  ©0  sin  ©.  »  (4-8) 

cos  Oj  sin  ©p  sin  ©Q)  -  cos  ©p  sin  ©Q 
a  -A  cos  ©.  cos  ©  f  A.  sin©  -A  sin  ©.  cos  © 

o  o  i  p  l  p  p  i  p 

Let  us  specify  that  the  magnitudes  of  the  misalignment  angles  {©j,  ©p,  ©0) 
be  limited  to  0.01  rad  (34.4'),  i.e.  |©.|  <0.01  rad,  where  j  =  i ,  p,  o.  Recall 

J  f 

that  these  misalignment  angles  Include  misalignments  of  the  mounting  fixture  and  of 
the  centrifuge  mounting  surfaces  with  respect  to  the  centripetal  and  gravity  vectors 
but  not  misalignments  due  to  angular  compliance.  We  may  now  use  the  small  angle 
approximations: 

-7 

sin  0^  ~  0^  with  an  absolute  error  less  than  3.3  x  10 

cc  .  o.  ^  1  -  4  o?  with  an  absolute  error  less  than  4  x  10  ^ 

J  2  j 

Substitute  the  above  approximations  in  Eqs.  (4-8)  and  neglect  third  order  errors, 
i.e.  triple  products  of  ©^,  j  =  i,  p,  o  and  products  of  ©^  and  the  compliance 
coefficients  k^,  k^,  kz> 

»i  *  A|«-5*p ‘  Ap«V*iV'Ao  Wi* 

aP  -  v'-K  -5*i> *  Ao'0it*pV-Ai#o  •  (4-9; 

aQ  =  A0  (1 -|©? -~©p)  +  Ai^p“Ap®i 
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SECTION  V 


MOUNTING  POSITIONS  AND  LIMITS 


In  order  to  determine  the  coefficients  of  the  model  equation,  it  is  necessary  to 
test  the  accelerometer  in  a  number  of  different  mounting  positions  and  at  many 
acceleration  levels.  The  chosen  positions  (others  are  possible)  are  shown  in  Figs.  3 
and  4.  Recall  that  the  centripetal  acceleration  vector  aQ  is  in  the  minus  x  direction 
and  the  effect  of  gravity  is  equivalent  to  an  upward  acceleration  of  one 

Table  I  shows  the  model  equation  coefficients  that  will  be  determined  in  each 
position  pair  of  Figs.  3  and  4. 


Table  I 


Position  Pair 

Coefficient 

1  and  2 

K2i! 

KC  K3i 

3  and  4 

V 

K3P 

5  and  6 

K2o’ 

K3o 

7  and  8 

Kip 

9  and  10 

Kpo 

11  and  12 

Koi 

Based  on  past  experience  with  this  type  of  torque  balanced  pendulous 
accelerometer,  we  may  conservatively  assume  that  the  model  equation  coefficients 
satisfy  the  following  limits : 


IK* I  •  I K2p !  ■  t K20 1  •  IV  ,  |K0, 1  ,  |K,|  <10"VS2 

lK9i!-IV-!K3ol  i10"4^ 

lkx  I  »  Iky  I  ,  !kz  I  <  10'°  rad/g  (1  arc  sec/g) 


(5-1) 
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(c)  POSITION  9 


(d)  POSITION  10 


Figure  4.  Accelerometer  Mounting  Positions  (Top  Views) 


Recall  that  we  have  already  limited  the  misalignment  angles  to 
|0j  |  <0.01  rad  (34.4*),  j  =  i,  p.  o. 

Further,  in  our  development  of  the  centrifuge  equations,  let  us  neglect 
errors  less  than  the  following  magnitudes: 

-0 

Bias  terms  <10  g 

— 

Linear  terms  <10  g/g 

-7  2 

Quadratic  terms  <  10  g/g 

o  3 

Cubic  terms  <  10  g/g 
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SECTION  VI 


INPUT  AXIS  NONLINEAR  COEFFICIENTS  K2i>  Kj., 


6.1  MOUNTING  AND  DATA  TAKING 

For  maximum  accuracy  in  the  separation  of  the  quadratic  and  cubic  terms, 
data  should  be  taken  with  the  input  acceleration  along  both  the  positive  and  the 
negative  input  reference  axis  such  as  positions  1  and  2  of  Fig.  3. 

Mount  the  accelerometer  on  its  fixture  in  position  1  and  check  for  proper 
operation.  Align  the  input  axis  (IA  not  IRA)  to  the  radial  direction  on  the  centrifuge 
within  0.01  rad  (34  min.)  about  the  pendulous  and  output  axes.  Align  the  output  axis 
(OA)  to  the  vertical  within  0.01  rad.  Determine  the  radius  arm  to  the  centroid  of 
the  proof  mass  by  whatever  procedure  is  appropriate  for  the  particular  facility. 

Measure  and  record  the  nominal  bias  and  scale  factor  by  the  two-point  static 
test  method  (IA  up  and  IA  down)  before  and  after  the  group  of  runs  in  each  position. 
Preferably,  these  should  be  made  in  the  centrifuge  environment.  The  measurements 
are  to  be  used  as  a  check  on  the  drift  or  shift  in  the  accelerometer  bias  and  scale 
factor  due  to  centrifuge  action  or  other  environmental  changes.  Any  drift  or  shift 
of  these  parameters  will  affect  the  accuracy  of  the  centrifuge  tests. 

Tne  accelerometer  should  be  tested  at  a  number  of  acceleration  levels  which 
are  approximately  evenly  spaced  over  the  input  range  in  each  mounting  position.  The 
centrifuge  speed  should  be  smoothly  and  unidirectionally  increased  from  one  acceleration 
level  to  the  next  and  allowed  to  stabilize  at  each  level  since  there  may  be  heating  effects 
at  the  different  torquing  levels.  When  the  accelerometer  has  been  tested  at  its  peak 
acceleration  level,  decrease  the  centrifuge  speed  smoothly  and  unidirectionally  from 
one  acceleration  level  to  the  next.  Test  at  approximately  the  same  levels  as  for 
increasing  speed.  In  each  mounting  position  the  accelerometer  should  be  tested  at 
a  minimum  of  11  acceleration  levels:  at  5  acceleration  levels  each  for  increasing 
and  decreasing  speeds  plus  the  peak  acceleration  level.  More  acceleration  levels  are 
desirable  since  redundancy  reduces  the  uncertainties  in  the  coefficients  determined  by 
the  centrifu  e  tests. 
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The  output  of  the  accelerometer  should  be  averaged  over  a  number  of 
revolutions  sufficient  to  minimize  quantization  errors;  quantization  may  be  a  function 
of  the  instrumentation  as  well  as  of  the  accelerometer.  The  period  of  the  centrifuge 
relative  to  earth,  the  output  of  the  accelerometer,  the  radius  arm,  the  arm  stretch, 
and  the  angular  compliance  (if  available)  must  be  determined  as  precisely  as  possible 
and  recorded.  It  is  common  practice  to  make  three  sets  of  measurements  at  each 
acceleration  level  and  to  use  either  the  median  or  the  mean  of  the  three  sets  of 
measurements  in  the  data  reduction. 

6.2  OUTPUT  EQUATIONS 

In  position  1  of  Fig.  3,  the  acceleration  components  along  the  input,  pendulous, 
and  output  reference  axes,  see  Eqs.  (4-7),  are: 


Substitute  the  above  values  in  Eqs.  (4-9).  Let  a  subscript  1  denote  the 
misalignments  and  outputs  in  position  1  and  a  subscript  r  denote  the  variable 
acceleration  level  for  this  position.  Let  n  equal  the  number  of  acceleration  test 
levels  in  position  1  so  that  r  =  1,  2,  3,  . . . ,  n.  From  the  limits  set  in  Eqs.  (5-1) 
and  (5-2),  it  is  evident  that  we  may  neglect  products  of  the  misalignment  angles  by 
the  compliance  coefficients,  e.g.  k^ 


ail  acr  (1  “  2  0  pi  "  2  0  ol  +  V  “  %!  “  0ol  0il) 


apl  acr  acr  <0ol  +  ^  +  ^il  +  ^pl  ^ol^ 

ao!  ~  acr  ky  +  acr  %1  +  ^  1  ~  2  ^il  ~  2  %1  ^ 


Substitute  Eqs.  (6-1)  in  the  model  equation,  Eq.  (2-1),  and  apply  the  limits  of 
Eqs.  (5-1)  and  (5-2).  The  acceleration  indicated  by  the  accelerometer  in  position  1 
is  given  by  Eq.  (6-2). 
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Alr  °  [K0  -  V  T  »ol  °il  +  K2o  +  K3d  +  Kpo  #il  '  Kot  “pi] 


+  a. 


1  .2 


cr 


[‘  -Ispl  -5  "Si  1  kv  -  <2K2i  +  2K.  -  2K2o  -  »V*pl 


+  Kip  *11  -  V  »ol  "  Koi  ]  *  “cr  [K2i  *  Kt  +  «Koi  -  3K3i>  “pi 
-  Kip  *ol]  *  K3i  “or 

In  position  2  of  Fig.  2  the  acceleration  components  along  the  accelerometer 
reference  axes  are: 

a12  M  “u  '  -  (1  *  V  “cs 
Ap2  '  av  =  kx  acs  -  kz  acs 


(6-2) 


Ao2  aw  1  ~  ky  acs 

Substitute  the  above  values  in  Eqs.  (4-9).  Let  a  subscript  2  denote  the 
misalignments  and  outputs  in  position  2  and  a  subscript  s  denote  the  variable 
acceleration  for  this  position.  Let  m  equal  the  number  of  acceleration  test  levels 
in  position  2  so  that  s  =  1,  2,  3,  . . . ,  m.  Note  that  the  misalignment  angles  with 
the  accelerometer  mounted  in  position  2  will  not  be  the  same  as  the  misalignment 
angles  with  the  accelerometer  mounted  in  position  1,  in  general. 


ai2  ~  "  acs  (  1  ~  2  0p2  "  2  0o2  +  ky  *  "  (0p2  "  0o2  0i2* 


ap2  *  kzV  +  acs  <0o2  +kx}  +  (0i2  +  0p2  0o2> 


o2 


“  acs  ky  “  acs  0p2  +  ^  1  “  2  0i2  ~  2  °p2^ 


(6-3) 


Substitute  Eqs.  (6-3)  in  Eq.  (2-1)  and  apply  the  limits  of  Eqs.  (5-1)  and  (5-2). 
The  indicated  acceleration  is  given  by  Eq.  (6-4). 
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(6-4) 


2s 


[X0  ■  °p2  *  °o2  °i2  +  K2o  *  K3o  +  Kpo  '>i2  '  Koi  0pa] 

'  “cs  t1  -  I  °p2  -  5  °o2  +  ky  -  <2K2i  *  2Kt  ‘  2K20  ’  31So>  *pS 


+  ^ip  ^i2  ^po°o2  ^oi  +  acs 


y)+  4  [K2i  ■  Kt 


+  <Koi  -  3K3i>  *p2 


“  Kip  0O2  ]  "  K3i 


cs 


It  is  evident  that  the  unknown  misalignment  angles  and  the  uncertainty  in  the 
radius  arm,  which  affects  the  value  of  a  and  a  ,  makes  it  impossible  to 

Ll  V/S 

precisely  determine  either  the  bias  or  the  scale  factor  from  Eq.  (6-2)  and/or 
Eq.  (6-4).  If  (K9i  +  Kt)  and  (K2i  -  Kfc)  are  to  be  determined  from  Eqs.  (6-2) 
and  (6-4)  with  an  errcr  less  than  10  percent  (aside  from  noise  in  the  data),  then 
the  magnitudes  of  the  terms  3K3i,  Kjp,  and  KQ.  may  not  exceed  the  magnitude 
°f  Kgj  or  Kj,  whichever  is  larger,  by  more  than  a  factor  of  two  or  three.  If  the 
magnitude  of  these  coefficients  exceed  those  limits  or  if  the  allowable  error  is  less 
than  10  percent,  then  the  misalignment  angles  (j  =  i,  o,  p)  must  be  reduced 
proportionately.  If  the  above  conditions  are  met,  then  Eqs.  (6-2)  and  (6-4)  may 
be  simplified  to: 


Alr  =  B0  +  B1  acr  +  B2  acr  +  K3i  acr 


A2s  C0  "  C1  acs  +  C2  acs  "  K3i  aes 


(6-5) 

(6-6) 


where 

B0,  CQ  =  constant  terms  of  Eqs.  (6-2)  and  (6-4),  respectively 
Bj,  Cj  =  linear  terms  of  Eqs.  (6-2)  and  (6-4),  respectively 

B2  ^  K2i  +  ^ 

C2  =  K2i  "  Kt 
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Then 


K2i  =  <B2  +  C2)/2 


(6-7) 


K(  =  (B2  -  C2)/2 


(6-8) 


The  effect  of  arm  stretch  *and  earth  rate  must  be  considered  in  determining 

a  and  a  ,  if  necessary.  Remember  that  both  a__  and  a  a.’e  the  centripetal 
cr  cs  cr  cs 

accelerations  and  are  always  positive. 

6.3  DETERMINATION  OF  THE  COEFFICIENTS 


In  the  proposed  data  reduction  procedure,  the  method  of  least  squares  is 
applied  to  all  the  data  at  once  rather  than  in  the  three  stage  method  of  Evans  and 
Fuhrman.  In  order  to  accurately  determine  the  coefficients  for  the  quadratic  and 
cubic  terms,  it  is  vital  to  have  the  cubic  coefficient,  at  least,  be  common  to  both 
positions  as  in  Eqs.  (6-5)  and  (6-6).  Let  the  best  fit  combined  outputs  of  positions  1 
and  2  as  obtained  from  Eqs.  (6-5)  and  (6-6)  be 


Ajk  =  <B0  +  Blacr  *  B2acr)6r  +  (C0  "  C1  acs  +  C2acs)Ss 


+  K3i<acr5r-acSV 


(6-9) 


where 

=  best  fit  value  of  output  for  jk  =  lr  or  2s 
Ajk  =  measured  value  of  output  for  jk  =  lr  or  2s 
6r  =  1  and  <5g  =  0  when  k  =  r. 

6g  =  1  and  6r  =  0  when  k  =  s. 


The  residual  at  centripetal  acceleration  level  afik  (k  =  r  or  s)  is; 
rk  =  Ajk  ‘  Ajk  =  (Alr  6r  +  A2s  6s)  "  [(B0  +  B1  acr  +  B2  acrJ  5r 

+  <C0  '  C1  acs  +  C2  acs>6s  +  K3i  <4  6r  “  acs  5s  >] 
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and  the  sum  of  the  squares  of  all  the  residuals  is 


(6-11) 


In  the  method  of  least  squares  (Ref.  3),  the  sum  of  the  squares  of  the 
residuals  is  minimized  by  setting  the  partial  derivatives  of  Eq.  (6-11)  with  respect 
to  each  unknown  coefficient  (Bq,  B^,  B9,  Cq,  C^,  C2,  and  K3i)  equal  to  zero.  Hie 
resulting  normal  equations  are  given  in  matrix  form  by  Eq.  (6-12).  Recall  that 
a  and  a  are  centripetal  accelerations  and  are  always  considered  to  be  positive. 

1  L  b 

Solve  Eqs.  (6-12)  and  (6-10)  for  the  unknown  coefficients  and  each  residual, 
respectively.  The  coefficients  BQ  and  CQ  should  closely  approximate  the  bias 
determined  from  a  two  or  four  position  test,  within  approximately  O.Olg.  The 
coefficients  B1  and  should  closely  approximate  unity  with  an  error  less  than 
0. 15  percent  if  the  radius  arm  has  been  determined  within  0. 1  percent.  The 
coefficients  K2i  and  Kj.  of  the  model  equation  are  determined  from  Eqs.  (6-7) 
and  (6-8).  A  plot  of  the  residuals  vs.  g  level  will  be  used  to  estimate  the  adequacy 
of  the  model  equation.  The  residuals  will  also  be  used  to  determine  the  uncertainty 
of  the  coefficients. 


If  =  0,  as  may  occur  in  puise  torqued  systems  of  constant  power  input, 


then  B2  =  C2  =  K2i  and  the  best  fit  equation  becomes: 


Ajk  ■  <B0  +  B1  acr>6r  +  <C0  '  °1  acS>6S  +  K2i<acr6r  +  acS6S» 


+  K3i<48r  '  acs6s> 


(6-13) 


24 


The  residuals  are: 


rk  “  Ajk'Ajk  =  (Alr6r*A2sds>- 


+  K^.ia"  6  +  a2  6  )  +  K„.{a3  6  -  a3  6  ) 
2iv  cr  r  cs  s'  3r  cr  r  cs  s' 


(B0  +  Blacr>6r  +  (C0-Ulacs)6s 

(6-14) 


The  normal  equations  are  obtained,  as  before,  by  taking  partial  derivatives 
of  the  sum  of  the  squares  of  the  residuals  with  respect  to  each  unknown  coefficient 
(BQ,  Bj,  CQ,  Cj^,  Kg-,  and  K^.)  and  setting  each  derivative  equal  to  zero.  The 
normal  equations  are  given  in  matrix  form  by  Eq.  (6-15). 

Solve  Eqs.  (6-15)  and  (6-14)  for  the  unknown  coefficients  and  each  residual, 
respectively. 

6.4  STANDARD  DEVIATION  OF  THE  RESIDUALS 

The  standard  deviation  of  the  residuals  is  a  measure  of  the  fit  of  the  mode’ 
equation  to  the  measured  accelerometer  output.  It  is  not  necessarily  a  measure  of 
the  accuracy  or  completeness  of  the  model  equation  in  describing  the  phenomena. 

The  unbiased  estimate  of  the  standard  deviation  of  the  residuals  is  given 
by  the  formula  of  Eq.  (6-16)  (Ref.  3). 


where 


$(*«,)  = 


/\ 

= 


z4  = 

n  +  m  = 

q 


k 


(n+m)  -  q 


unbiased  estimate  of  the  standard  deviation  of  the 
residuals  r^  (k  =  r,  s)  -  £ 

o 

sum  of  the  squares  of  all  the  residuals  -  g 

total  number  of  data  points  for  k  equal  r  and  s 

total  number  of  coefficients  Bq,  Cq,  etc. 

E.g.,  q  =  7  in  Eq.  (6-12) 


(6-16) 
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6.5 


UNCERTAINTY  OF  THE  COEFFICIENTS 


The  unbiased  estimate  of  the  uncertainty  of  the  coefficients  (Ref.  3)  obtained 
from  Eq.  (6-12)  or  (6-15)  are 

£(Q>  ^  .(^  (6-17) 

where 

e(Q)  -  uncertainty  of  coefficient  Q 

Q  -  one  of  the  coefficients  in  the  second  (column)  matrix  of 

Eq.  (6-12)  or  (6-15),  such  as  Bg 

D  -  determinant  of  the  left-hand  matrix  in  Eq.  (6-12)  or  (6-15) 

M  (Q)  -  minor  of  the  determinant  D  for  the  coefficient  Q.  E.g. , 

M(B2)  is  obtained  by  deleting  column  and  row  indicated  by 
asterisks  in  Eq.  (6-12) 

It  is  quite  possible  that  a  coefficient  of  Eq.  (6-12),  or  of  a  similar  equation, 
should  actually  be  zero  but,  because  of  the  small  number  of  observations,  it  is 
practically  certain  that  we  would  obtain  a  nonzero  value.  We  may  determine  if  a 
coefficient  Q  is  likely  to  be  zero  by  assuming  the  null  hypothesis,  i.e. ,  Q  =  0, 
and  applying  a  statistical  significance  test,  such  as  the  Student’s  t  statistic  to  the 
significance  ratio  |Q|  /<r(Q).  The  Student’s  t  statistic  varies  with  the  number  of 
degrees  of  freedom  (f  =  o-un-q)  and  the  chosen  rejection  level  for  the  null  hypothesis. 
The  rejection  level  depends  on  the  application  and  is  a  matter  of  judgment;  however, 
a  5  percent  rejection  level  is  commonly  used.  By  a  5  percent  rejection  level,  we 
mean  that  the  null  hypothesis  will  be  rejected  if  there  is  only  a  5  percent  chance 
that  Q  =  0  if  |Q|  /  a  (Q)  of  our  sample  exceeds  a  certain  value.  This  is  equivalent 
to  saying  that  we  can  be  95  percent  confident  that  the  coefficient  is  significant  if  the 
significance  ratio  exceeds  the  Student's  t  statistic. 

At  the  5  percent  rejection  i.evei  (2  sided),  the  t  statistic  varies  from  3.182 
for  f  -  3  to  1.96  for  f  -  co .  Since  we  should  have  a  minimum  of  22  data  points 
(n+m  2>  22)  and  since  the  maximum  number  of  coefficients  ir  Eq.  (6-12)  is  7  (q  =  7), 
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the  minimum  number  of  degrees  of  freedom  is  f  =  n+m-q  -  15.  At  the  5  percent 
rejection  level  and  f  =  15,  the  Student’s  t  statistic  is  2.13.  With  our  assumptions, 
the  t  statistic  can  vary  only  from  2. 13  for  f  =  15  to  1.96  for  f  =  oo  .  Therefore,  we 
may  apply  the  simple  rule  that  the  null  nypothesis  will  be  rejected  if  |Q|  /a  (Q)  >  2. 
For  a  more  extensive  discussion  on  the  null  hypothesis  and  the  various  significance 
tests,  see  Ref.  4,  or  equivalent. 

If  a  coefficient  is  set  equal  to  zero  as  a  result  of  the  t  test,  Eqs.  (6-2)  and 
(6-4)  thru  (6-15)  should  be  modified  appropriately. 

6.6  "GOODNESS”  OF  THE  MODEL  EQUATION 

6.6.1  Introduction 

There  are  no  general  rules  for  defining  the  "goodness"  of  a  model  equation 
though  we  may  apply  one  or  more  criteria  such  as: 

(1)  Is  it  derived  from  basic  physical  principles? 

(2)  Is  it  simple? 

(3)  How  well  does  it  fit  the  observational  data  ? 

(a)  What  is  the  value  of  the  standard  deviation  of  the  residuals  ? 

(b)  What  is  the  magnitude  of  the  peak  residual  ? 

r' 

(4)  Would  a  different  model  equation  fit  the  observational  data  better? 

It  is  apparent  that  the  "goodness"  of  a  model  equation  is  a  subjective  question 
which  caii  be  settled  only  by  each  person  for  each  application.  There  are  some 
statistical  criteria  that  we  may  apply,  but  they  must  be  applied  with  judgment. 

TTius,  in  Section  6.5,  we  showed  how  to  apply  the  null  hypothesis  and  the 
Student’s  t  significance  test  to  the  coefficients  of  our  best-fit  equation,  but  at  an 
arbitrary  5  percent  rejection  level,  which  was  a  matter  of  judgment.  In  examples  1 
and  2,  below,  we  will  illustrate  some  of  the  decisions  or  judgments  that  may  be 
required  in  adopting  a  model  equation. 


29 


It  is  obvious  that  observational  data  must  bo  valid,  i.o.  without  systematic 
error  and  with  a  minimum  of  noise,  if  we  are  to  obtain  a  "good"  model  equation. 

Care  must  be  taken  in  the  choice  of  instrumentation,  the  test  setup  (including  temperature 
control),  the  mounting  fixture(s),  the  data-taking  procedures,  and  data  processing 
(such  as  smoothing)  in  order  to  avoid  systematic  errors  and  excessive  noise. 

Use  of  one  model  equation  for  all  accelerometers  in  all  applications  is  to 
be  deplored.  A  model  equation  which  is  satisfactory  for  a  particular  model  of 
accelerometer,  in  a  particular  application,  may  be  quite  unsatisfactory  for  a 
different  accelerometer  and/or  a  different  application.  Equally  deplorable  is  a  data 
reduction  program  which  grinds  out  coefficients,  even  programs  with  statistical 
criteria,  without  the  intervention  of  human  interpretation  and  judgment. 

6.6.2  Example  1 

The  effects  of  an  incorrect  model  equation  will  now  be  illustrated  with  an 
artificial  example.  For  acceleration  inputs  only  along  plus  and  minus  IA,  let  the 
true  and  exact  output  equation  be: 

A  -  2xl0“3  +  a.  +  10‘°  af  +  10-fi  a?  +  lO-7  a?  (6-18) 

l  i  t  l  '  7 

The  calculated  outputs  A  over  the  input  range  from  +30  £  to  -30  g  in  2  £ 
steps  is  given  in  column  2  of  Table  II  with  the  results  rounded  off  at  six  decimal 
places.  'Hie  tabulated  inputs  and  outputs  of  columns  1  and  2  will  now  be  treated  as 
observational  data  to  which  we  will  fit  a  cubic  equation  by  the  method  of  least  squares. 

Ihe  resulting  cubic  equation  is: 

A*  -  -5.83326  x  10~3  +  a.  +  9.2xl0-5  a3  +  10"6  a?  (6-19) 

Note  that  the  linear  and  cubic  coefficients  are  unchanged,  as  would  be 
expected,  since  it  is  only  the  symmetric  fourth  degree  term  which  has  been  omitted. 

Also  note  that  the  coefficient  is  an  order  of  magnitude  larger  than  that  of 
Eq.  (6-18). 
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Table  II 


1 

2 

3 

4 

5 

a. 

1 

A 

A' 

r'  =  A’ -A 

i< 

;i 

~u 

30 

30. 119000 

30. 10396674 

-0.01503326 

-0.007200 

28 

28.093258 

28.08824674 

-0. 00501126 

0.002822 

26 

26.072034 

26.07393474 

0.00190074 

0.009734 

24 

24.054762 

24. 06098274 

0.00622074 

0.014054 

22 

22.040914 

22.04934274 

0. 00842874 

0.016262 

20 

20.030000 

20. 03896674 

0. 00896674 

0.016800 

18 

18.021570 

18.02980674 

0.00823674 

0.016070 

16 

16.015210 

16.02181474 

0. 00660474 

0.014438 

14 

14.010546 

14. 01494274 

0.00439674 

0.012230 

12 

12.0072*2 

12. 00914274 

0. 00190074 

0,009734 

10 

10.005000 

10.00436674 

-0. 00063326 

0.007200 

8 

8.003562 

8.00056674 

-0.00299526 

0.004838 

6 

6.0Q2706 

5.99769474 

-0.00501126 

0.002822 

4 

4.002250 

3.99570274 

-0.00654726 

0.001286 

2 

2. 002050 

1.99454274 

-0.00750726 

0.000326 

0 

0.002000 

-0.00583326 

-0. 00783326 

-2 

-1.997966 

-2.00547326 

-0. 00750726 

0.000326 

-4 

-3.997878 

-4,00442526 

-0.00654726 

0.001286 

-6 

-5.997726 

-6.00273726 

-0.00501126 

0.002822 

-8 

-7.997462 

-8.00045726 

-0.00299526 

0.004838 

-10 

-9.997000 

-9.99763326 

-0.00063326 

0.007200 

-12 

-11.996214 

-11.99431326 

0. 00190074 

0.009734 

-14 

-13.994942 

-13.99054526 

0.00439674 

0.012230 

-16 

-15.992982 

-15.98637726 

0. 00660474 

0.014438 

-18 

-17.990094 

-17.98185726 

0.00823674 

0.016070 

-20 

-19.986000 

-19.97703326 

0.00896674 

0.016800 

-22 

-21.980382 

-21.97195326 

0.00842874 

0.016262 

-24 

-23.972886 

-23.96666526 

0. 00622074 

0.014054 

-26 

-25.963118 

-25.96121726 

0. 00190074 

0.009734 

-28 

-27.950646 

-27.95565726 

-0.00501126 

0.002822 

-30 

-29.935000 

-29.95003326 

-0.01503326 

-0.007200 

Rms  Value  = 

6.9  x  10"3 

10. 5  x  lo"3 

£r'  =  lxlO"6 

V 

r 

s' 

I 

l 


b 

y 


The  outputs,  A\  and  the  corresponding  residuals,  A*  -  A.  were  calculated 
to  8  decimal  places  and  are  tabulated  in  columns  3  and  4  of  Table  Ii.  The  sum  of  the 
residuals  is  not  quite  zero  due  to  round-off  errors. 

Using  Eq.  (6-1(5),  the  unbiased  estimate  of  the  standard  deviation  of  the 
residuals  is  a  (r*)  -  7.422  x  10-3  g.  The  coefficients  and  the  uncertainties  of  the 
coefficients  of  Eq.  (6-19)  are  listed  along  with  the  significance  ratio  |Q|  /  a  (Q)  in 
Table  III. 

Table  III 


Q 

a  (Q) 

i  Q 1  /  (Q) 

KQ  =  -5. 83326xl0"3 

1.333xl0~3 

4.37 

K'i  *  1 

0.1S7xl0“3 

5.35 

K9  -  9.2xl0_b 

0.467x10“° 

19.7 

k3  -  ixuf6 

0.298xl0_fi 

3.36 

Note  that  K'j  is  not  the  scale  factor  as  defined  in  Eq.  (2-1)  but  is  in  g/g. 

The  application  of  the  Student's  t  statistic  for  f  =  31-4=27  degrees  of  freedom 
shows  that  all  four  coefficients  are  significantly  different  than  zero  at  a  99  percent  or 
higher  confidence  level.  (1  percent  rejection  level)  However,  two  of  the  coefficients 
are  far  different  than  the  corresponding  coefficients  of  the  exact  model  equation, 

Eq.  (6-18). 

Let  us  not  stop  here  but  let  us  plot  the  residuals  given  in  the  fourth  column 
of  Tabic  II,  sec  Fig.  5.  Even  a  casual  glance  shows  that  the  residuals  are  not 
randomly  distributed,  therefore,  the  residuals  are  not  primarily  due  to  noise  in  the 
observational  data  but  arc  largely  systematic.  From  the  shape  of  the  curve  (3  maxima), 
we  should  suspect  that  a  fourth  degree  term  should  be  added  to  the  model  equation. 

If  this  is  done,  wt*  know  that  we  would  get  Eq.  (6-18)  except  for  slight  computational 
errors.  If  there  had  been  four  maxima  in  Fig.  5  instead  of  three,  we  would  have 
suspected  that  the  model  equation  required  a  fifth  degree  term. 
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Figure  5.  Residuals,  Ex.  1 
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It  should  be  remembered  that  this  illustration  is  an  idea)  case  with  no  noise 
in  the  data  except  that  due  to  rounding-off  errors.  In  an  actual  case,  the  residual  plot 
would  be  distorted  and  roughened  by  noise. 

Let  us  now  consider  one  more  fact  that  is  best  illustrated  by  this  artificial 
example.  As  indicated  previously,  it  is  presumed  that  the  bias  and  scale  factor  obtained 
by  reduction  of  centrifuge  data  is  probably  seriously  in  error  due  to  rr  isalignments 
and  uncertainty  in  the  radius  arm  to  the  centroid  of  the  proof  mass.  The  bias  and 
scale  factor  are  generally  obtained  from  a  two,  four,  or  six  point  static  calibration 
test  in  a  one  £  field.  In  this  example  we  would  find  it  to  be  those  of  Eq,  (6-18),  i.e. 

K0  =  0. 002 £  and  K'j  =  1  g/j£.  The  nonlinear  coefficients  would  be  those  of  Eq.  (6-19). 
Thus,  it  would  be  presumed  that  the  output  equation  is: 

A"  =  2  x  10"3  +  a.  +  9.2xl0“5  a?  +  l(f6  a3  (6-20) 

In  many  applications,  the  bias  is  subtracted  from  the  output  but  the  other 
errors  are  tolerated;  the  errors,  for  this  example,  would  be: 

r  =  A"  -  (a.  +  0.002)  (6-21) 

In  other  applications,  not  only  the  bias  but  also  the  errors  due  to  the 
nonlinear  terms  are  removed  by  a  computer.  The  errors  would  be  the  difference 
between  the  predicted  output  of  Eq.  (6-20)  and  the  true  output  of  Eq.  (6-18),  i.e. 
the  errors  would  be: 

r"  =  A"  -  A  (6-22) 

The  residuals  r"  are  given  in  column  5  of  Table  II.  As  would  be  expected, 
the  rms  value  of  r"  is  larger  than  the  rms  value  of  r':  in  this  example,  over 
50  percent  larger. 

Example  1  illustrates  that  it  is  not  sufficient  to  look  at  the  standard  deviation 
of  the  residuals  nor  to  compute  the  statistical  significance  of  the  coefficients  to 
determine  if  we  have  a  "good"  model  equation;  we  must  also  look  at  a  plot  of  the 
residuals.  If  the  residuals  look  as  if  they  are  randomly  distributed,  we  probably 
cannot  obtain  a  better  fit  model  equation,  but  we  might  try  a  polynomial  of  lower 
degree  to  determine  if  it  is  satisfactory.  If  the  residuals  are  not  randomly  distributed, 
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we  might  try  a  more  complex  model  equation,  or  we  might  decide  that  the  model 
equation  is  satisfactory  for  the  particular  application,  i.e.  the  errors  are  acceptable. 
However,  we  should  also  recognize  that  the  actual  residuals  are  larger  than  the 
calculated  residuals  as  illustrated  by  the  difference  between  columns  5  and  4  of 
Table  II. 

6.6.3  Example  2 

In  this  example,  we  will  use  actual  experimental  data  and  obtain  certain 
coefficients  of  the  model  equation,  Eq.  (2-1):  first,  with  the  assumption  Kj.  i  0 
and  second  with  the  assumption  that  Kt  =  0.  The  input  and  output  data  of  the  binary 
pulse  torque  accelerometer  (constant  torquing  power)  are  listed  in  Tables  IV  and  V. 
The  input  data  have  been  partly  processed  to  save  time  and  space,  i.e.  the  input 
acceleration  has  been  derived  from  the  period  of  the  centrifuge  and  the  radius  arm 
with  corrections  for  arm  stretch  and  earth  rate.  At  the  end  of  the  increasing  run, 
the  centrifuge  was  run  beyond  45  g  before  starting  the  decreasing  run. 

The  data  in  Tables  IV  and  V  are  substituted  in  Eqs.  (6-12)  and  (6-15)  to 
determine  the  coefficients  BQ,  CQ,  etc.  The  resulting  coefficients,  the  unbiased 
estimate  of  the  standard  deviation  of  the  residuals  from  Eq.  (6-16),  the  uncertainties 
of  the  coefficients  from  Eq.  (6-17),  the  number  of  degrees  of  freedom,  and  the 
significance  ratios  are  given  in  Table  VI.  The  residuals  were  determined  using 
Eqs.  (6-10)  and  (6-14)  and  are  plotted  in  Figs.  6  and  7. 

Normally,  we  would  expect  the  unbiased  estimate  of  the  standard  deviation  of 
the  residuals  for  Eq.  (6-12)  to  be  less  than  for  Eq.  (6-15)  since  Eq.  (6-12)  has  more 
coefficients  and,  therefore,  it  should  give  a  better  fit.  Actually,  the  rms  of  the 
residuals  is  less  for  Eq.  (6-12)  by  approximately  1  percent  but  this  is  more  than 
offset  by  the  effect  of  the  reduction  in  the  degrees  of  freedom  from  30  to  29. 

'Ihe  significance  ratio  for  B2>  Cg,  and  Kgi  of  Eq.  (6-12)  indicate  that  these 
coefficients  may  not  be  significantly  different  from  zero.  On  the  other  hand,  all  the 
coefficients  of  Eq.  (6-15)  are  significantly  different  from  zero  at  a  99  percent  or 
higher  confidence  level.  This  is  not  too  surprising  when  we  consider  that  B2  and  C2 
are  of  opposite  sign  and  of  different  magnitudes;  therefore,  the  combination  has  some 
of  the  characteristics  of  both  quadratic  and  cubic  terms.  Thus,  we  have  three  terms 


Table  IV 


POSITION  1  OF  FIGURE  3  (at  >  0) 


Increasing  Data 

Decreasing  Data 

Input 

Output 

Input 

Output 

a 

A, 

a 

A. 

cr 

lr 

cr 

lr 

5.000281 

5.034095 

45.048050 

45. 119699 

10.020010 

10.059407 

40.040700 

40.108608 

15.018229 

15.063017 

35.022121 

35.086046 

20.030342 

20.080320 

30.020392 

30.080080 

25.015052 

25.069831 

25.015060 

25. 070161 

30.020392 

30.079824 

20.030342 

20.080584 

35.022121 

35.085875 

15.018229 

15.063135 

40.040700 

40.105535 

10.020008 

10.059338 

45.048050 

45.119553 

5.000280 

5. 033769 

Tabic  V 


POSITION  2  OF  FIGURE  3  (aj  <  0) 


1  . .  "  1  1  1  1  '  ■■  1 

Increasing  Data 

Decreasing  Data 

Input 

Output 

Input 

Output 

a 

A„ 

a 

A„ 

cs 

2s 

cs 

2s 

5.000281 

-4.976889 

45.048050 

-45.055412 

10.020010 

-10.000542 

40.040717 

-40.044992 

15.018225 

-15.002753 

35.022108 

-35.022991 

20.030336 

-20.018835 

30.020392 

-30.018014 

25.015060 

-25.007460 

25.015052 

-25.008900 

30.020392 

-30.016623 

20.030336 

-20.C20096 

35.022108 

-35.022125 

15.018229 

-15.003737 

40.040717 

-40.044235 

10.020010 

-10.001131 

45.048050 

-45.054973 

5.000281 

\ 

-4.976946 
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performing  the  function  of  two  terms,  at  least  partially,  and  all  three  are  weakened 
by  this  process.  In  other  words,  it  is  difficult  to  separate  an  odd  square  term  (a 
quadratic  term  with  change  of  sign  at  0)  such  as  the  | a .  |  a.  term  from  a  cubic  term. 

It  is  apparent  that  Eq.  (6-15)  is  a  more  appropriate  equation  to  use  than 
Eq.  (6-12)  for  this  example  since  it  is  simpler,  the  unbiased  estimate  of  the  standard 
deviation  of  the  residuals  is  smaller  (though  not  significantly  so),  and  the  coefficients 
are  statistically  more  significant  as  shown  by  the  significance  ratios. 

Though  we  have  shown  that  Eq.  (6-15)  is  better  than  Eq.  (6-12)  for  this  case, 
we  have  not  demonstrated  that  it  is  the  "best”  equation.  The  plots  of  the  residuals 
in  Figs.  6  and  7  confirm  our  previous  conclusion  that  Eqs.  (6-12)  and  (6-15)  fit  the 
data  about  equally  well  and  they  also  show  that  the  residuals  are  definitely  not  random. 

However,  unlike  our  Example  1,  we  cannot  reduce  the  residuals  significantly 
by  going  to  a  higher-order  polynomial  since  the  residual  errors  are  primarily  due  to 
hysteresis  for  negative  input  accelerations. 

At  this  time,  the  mechanism  causing  the  hysteresis  is  unknown.  The  hysteresis 
could  be  due  to  a  number  of  causes  of  which  some  may  be  inherent  to  the  accelerometer 
design  or  it  could  be  due  to  the  fixturing  (including  the  centrifuge)  or  to  the  instrumenta¬ 
tion. 


If  we  had  blindly  used  Eq.  (6-12)  and  accepted  the  derived  coefficients,  we 
would  have  determined  values  of  and  by  means  of  Eqs.  (6-7)  and  ("'-8). 

These  would  be  -1.412  ix  g/g  and  -2.76.1  fx  g/g2,  respectively.  However,  when  we 
applied  the  Student's  t  statistic  to  the  significance  ratio,  we  showed  that,  statistically, 
those  coefficients  were  not  significantly  different  than  zero.  Physically,  we  should 
not  be  surprised  if  K^.  =  0  since  this  is  a  binary,  pulse-torqued  accelerometer  with 
constant  power  input. 


**  *  lvnoisaa 


ZO  864SZ  U 
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INCREASING  DATA 
DECREASING  DATA 


Residuals  of  Example  2  with  Kt  =  0,  Eq.  (6-15) 


SECTION  VII 


PENDULOUS  AXIS,  NONLINEAR  COEFFICIENT 


AND  K 


3p 


Positions  3  and  4  of  Fig.  3  will  be  used  to  obtain  the  value  of  K2p  and  Kgp.  The 
mounting  and  data  taking  procedures  for  the  accelerometer  in  these  positions  would  be 
similar  to  those  outlined  in  Section  6. 1. 


In  position  3  of  Fig.  3,  the  acceleration  components  along  the  accelerometer 
reference  axes  are: 


a  +  k 
cr  z 


Ap3=-V<1  +  ky>acr 


A  _  =  a  =  1  -  k  a 
o3  w  y  cr 


Substitute  the  above  values  in  Eqs.  (4-9).  Let  a  subscript  3  denote  the  outputs 
in  position  3  and  a  subscript  r  denote  the  variable  acceleration  level  for  this  position. 
Let  n  equal  the  number  of  acceleration  test  levels  in  position  3  so  that  r  =  1,  2,  3, 

. . .  n.  Neglect  second  order  errors  as  was  done  in  Section  6. 2.  From  here  on,  it  is 
to  be  understood  that  the  misalignments  0.  (j  =  i,  p,  o)  will  be  different,  in  general, 
for  each  mounting  position.  For  convenience,  the  subscript  numbers  on  the  mis¬ 
alignments  will  be  dropped. 

ai3  *  acr  kz  f  acr  K  +  *1  »p  +  ‘  <*p  '  0o  ®i> 

ap3  ‘  acr  d  -  \  ®o  - \  A  +  V  +  (0i  +  W  '  <7-‘> 

ao3  =  -4ky.-acr''i  +  (1-5®f-5®p> 

J 

Substitute  Eqs.  (7-1)  in  the  model  equation  and  apply  the  limits  of  Eqfc’,  (5-1) 
and  (5-2).  Since  the  bias  and  linear  terms  are  so  uncertain  because  of  misalignments, 


Preceding  page  blank 


43 


etc. ,  we  will  not  list  the  individual  factors  in  these  terms.  The  acceleration  indicated 
by  the  accelerometer  may  be  expressed  in  the  following  form: 


A3r  "  D0  +  D1  "cr  + 


K2P  +  Kip  %  *  <3K3p  *  Kpo>°i  +  kz]  4  *  K3pacr  <7'2) 


In  position  4  of  Fig.  3,  the  acceleration  components  along  the  accelerometer 
reference  axes  are: 


A. .  -  -a  -  (k  a  +  k  a^  ) 
i4  v  *  x  cs  z  cs' 

p4  u  '  y  cs 

9 

A. -a  =  1  -  k  a“ 
o4  w  y  cs 

Substitute  the  above  values  in  Eqs.  (4-9).  Let  the  subscript  4  denote  the  outputs 
in  position  4  and  a  subscript  s  denote  the  variable  acceleration  level  for  this  position. 
Let  m  equal  the  number  of  acceleration  test  levels  in  position  4  so  that  s  =  1,  2,  3, 

. . . ,  m.  Neglect  second  order  errors.  The  acceleration  components  along  the  true 
axes  are: 


a. .  =  -a  k  -  a  (o  +  o.  o  +  k  )  -  (p  -  p  o.) 

[4  >7  no  '  a  ;  n  v'  r\ 


CS  Z  CS  x  O  ip 


owi' 


1  9  12 

a  ,  -  -a  _  (1  -  v  0  -  —  o.  +  k  )  +  (p.  +  0  0  ) 
p4  cs  x  2  o  2  l  y’  '  l  p  o’ 

a  .  =  -a2  k  +a  o.  +  (1  -  ^  o2  -  ^  02  ) 
o4  cs  y  cs  l  '  2  l  2  p  ’ 


(7-3) 


Substitute  Eqs.  (7-3)  in  Eq.  (2-1)  and  apply  the  limits  of  Eqs.  (5-1)  and  (5-2). 
The  indicated  acceleration  in  position  4  is 


A  .  =  E  -  E.  a  + 
4s  0  1  cs 


K2P  +  Kip  eo  +  <3K3P  "  V  »l  '  kz  ]  acs  -  K3p  4  <7'4> 


If  Kgp  is  to  be  obtained  from  Eqs.  (7-2)  and  (7-4)  with  an  error  of  no  more  than 
10  percent  (aside  from  noise  in  the  data),  then  the  magnitudes  of  the  products  of  the 
accelerometer  coefficients  by  the  misalignment  angles  which  appear  in  those  equations 
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If  those  products  are  sufficiently  small,  then  Eqs.  (7-2) 


must  be  smaller  than  K„  . 

2p 

and  (7-4)  may  be  simplified. 


A_  -  DA  +  D,  a  +  D„  a2  +  K0  a3 
3r  0  1  cr  2  cr  3p  cr 

A .  rE.-D.a  +E„a2  -  K0  a3 
4s  0  1  cs  2  cs  3p  cs 


where 


(7-5) 

(7-6) 


°2  =  K2P  +  k: 


The  coefficients  DQ,  EQ,  D.,  etc. ,  and  their  uncertainties  are  determined  by 
the  methods  outlined  in  Section  6.  If  it  is  assumed  that  the  compliance  coefficient  k 

z 

is  the  same  for  both  positions  3  and  4  (depends  on  construction  of  mounting  fixture) 
then 


K2p"f<D2  +  E2>  <7-7> 

The  residuals  should  be  determined  and  plotted  to  see  if  there  are  systematic 
errors  that  indicate  the  model  equation  is  incomplete.  The  standard  deviation  of  the 
residuals  and  the  uncertainty  of  the  coefficients  are  obtained  in  a  manner  similar  to 
that  outlined  in  Sections  6. 4  and  6, 5. 

Suppose  that  in  place  of  positions  3  and  4  of  Fig.  3,  we  were  to  use  the 
alternate  positions  of  3'  and  4'  of  Fig.  8  for  the  determination  of  and  K^.  The 
output  for  the  latter  positions  are: 


A3r=DitDlacr  + 


A’  =  E'  -  E,a  + 
4s  0  1  cs 


[  K2P  +  <K  ip  +  SV  °o  -  Kpo  0o  +  ky  ]  4  +  K3p4 

[  K2p  +  <Kip  -  3K3p>  %  '  Kpo  »i-ky]  4  -  K3p4 


(7-8) 

(7-9) 
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w,  IRA 


Figure  8.  Alternate  Mounting  Positions 


The  only  significant  difference  between  the  above  equations  and  Eqs.  (7-2)  and 
(7-4)  is  that  the  compliance  coefficient  k  has  replaced  the  coefficient  k  in  the 

y 

quadratic  term.  It  would  be  desirable  to  use  Eqs.  (7-8)  and  (7-9)  instead  of  Eqs.  (7-2) 
and  (7-4)  if: 


(a) 


k  is  measured  or  known  but 

y 


k 


z 


is  not 


(b) 

(c) 


if  k  <  k 

y  z 


if  k  is  the  same  for  positions  S’  and  4'  but  k  is  not  the  same  for 
y  z 

positions  3  and  4. 
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SECTION  VIII 


OUTPUT  AXIS,  NONLINEAR  COEFFICIENT  K2q  AND  K3o 


Positions  5  and  6  of  Fig.  3  will  be  used  to  obtain  the  value  of  K2q  and  K3q. 

The  mounting  and  data  taking  procedures  for  the  accelerometer  in  these  positions  would 
be  similar  to  thouse  outlined  in  Section  6. 1, 

In  position  5  of  Fig.  3,  the  acceleration  components  along  the  accelerometer 
reference  axes  are: 

A.c  =  a  =  1  -  k  a2 
i5  w  y  cr 

■  .  2 
ap5  av  "  xacr  +  z  acr 

Ao5  -  =  I1  +  ky>  acr 

substitute  the  above  values  in  Eqs.  (4-9).  Let  a  subscript  5  denote  the  outputs 
in  position  5  and  a  subscript  r  denote  the  variable  acceleration  level  for  this  position. 
Again  neglect  second  order  errors. 

9  19  12s 

A.„  =  -a  k  -  a  (<b  -  0  0,)  +  (l  -  4  0  -  t  0  )  I 
i5  cr  y  cr  ""p  vo  vi’  '  2  *p  2  vo '  f 


a  -  =  a  k  +  a  (0,  +  0  0  +  k  )  -  0 
p5  cr  z  cr  'vi  vp  vo  x;  vo 

t  12  12 

a  .  =  a  (1  "  t  0.  “70  +  k  )  +  0 
00  cr '  2  yi  2  vp  y'  p 


(8-1) 


Substitute  Eqs.  (8-1)  in  the  model  equation  and  apply  the  limits  of  Eqs.  (5-1) 
and  (5-2). 


A5r=F0+F 


lacr  +  I  K2o  +  Kpo  h  *  <3K3o  -  Kol>  V  kyl  &  +  K3o4  ^ 


In  order  to  minimize  the  uncertainty  in  K2o  introduced  by  the  compliance  k^, 
position  b  has  been  chosen  which  is  obtained  by  rotating  about  the  v  axis  rather  than 
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the  w  axis  as  was  done  to  obtain  positions  2  and  4.  The  acceleration  components 
aiop.g  the  accelerometer  reference  axes  are: 

A...  -a  =  -  (1  -  k  a2  ) 
i(>  w  '  y  cs 

o 

A  -  a  =  k  a  k  a" 
p‘o  v  x  cs  z  cs 

A  a  -  -  ( 1  +  k  )  a  „ 
o(>  u  y  cs 

In  the  usual  manner,  we  find  the  acceleration  components  along  the  true 
axes  are: 


o  ic>  12 

a.„  -  a”  k  +  a  (o  -  o  o.)  -  (1  -  —  0  -  —  0  ) 

lG  cs  y  cs  '  p  or'  2  p  2  0  ' 

V  =  i  k z  -  acs  <V  Vo  ‘  V  +  %  "  <8_3> 

ao(i  ■  -1cs  (1  -  2  °f  •  2  °p  *  V  ‘  °P 

Substitute  Eqs,  (8-3)  in  Eq.  (2-1)  and  apply  the  limits  of  Eqs.  (5-1)  and  (5-2), 

^6s  =  G0  ~  G1  acs  +  [K2o  +  KPo0i"<3K3o  +  Koi>l’p+kyj4-K3o4  <8'4> 

If  K2q  is  to  be  obtained  from  Eqs.  (8-2)  and  (8-4)  with  an  error  of  no  more 
than  10  percent  (aside  from  noise  in  the  data),  then  the  magnitudes  of  the  products  of 
the  accelerometer  coefficients  by  the  misalignment  angles  which  appear  in  those 
equations  must  be  much  smaller  than  K2q.  If  those  products  are  sufficiently  small, 
then  Eqs.  (8-2)  and  (8-4)  may  be  simplified. 


A5r“  F0  +  F1  acr  +  F2  4  +  K3o4 


AGs  ~  °0~°1  acs  +  °2  acs  ■  K3o  acs 


where 


F2  =  K2o  -  ky 


C2  *  K2o  +  ky 


■13 


Assuming  that 


is  the  same  for  positions  3  and  4,  then 


K 


2o 


= 


+  G2> 


(8-7) 


The  values  of  K„  and  K„„  are  obtained  in  a  manner  similar  to  that  used  for 
2o  3o 

the  determination  of  K„_  and  K„  in  Section  7, 

op 

If  there  is  any  advantage  to  having  the  compliance  coefficient  kz  instead  of 
in  the  quadratic  term  use  alternate  positions  in  which  IRA  is  along  the  plus  v  axis  when 
ORA  is  along  the  minus  u  axis  and  in  the  second  position  IRA  is  along  the  minus  v  axis 
when  ORA  is  along  the  positive  u  axis. 
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SECTION  IX 


CROSS  COUPLING  COEFFICIENT  K. 

ip 


For  the  assumed  model  equation,  the  best  mounting  positions  for  determining 

K.  are  those  in  which  the  acceleration  components  along  IRA  and  PRA  are  both 
»P 

reversed  in  sign.  This  results  in  a  common  cubic  term  which  is  vital  for  accuracy 

in  separating  the  quadratic  and  cubic  terms  of  the  output  equations.  In  order  to _ _ 

separate  the  compliance  term  /"i  kz  from  Kipf  it  is  necessary  to  choose  positions 
in  which  one  and  only  one  of  these  terms  reverse  sign.  Positions  7  and  8  of  Fig.  4 
is  one  set  of  positions  which  satisfy  the  above  criteria  providing  k  is  the  same  for 
the  two  positions. 

The  acceleration  components  along  the  accelerometer  reference  axes  in  these 
positions  are: 


“  r~  {an  *'  K)  i~  I  acr  kz  +  acr  ^  +  ^“x  +  ky) 


i7  yr*  v  yi" 


V  (au  +  V =  [aorkz-acr>1-kx  +  k,)J 


Ao7  aw  1  '  ky  acr 


Ai8  ~  y^  (au  ”  V 


=  ^[acskz  +  ara.(1+vv] 

Ap8  “  ^au  +  *V  ~  [acs  kz  "  acs  ^  “  w] 


ao4  aw  ~  1  ”  ky  acs 


(9-i, 


(9-2) 


Substitute  Eqs.  (9-l),-and  (9-2)  in  Eqs.  (4-9)  to  obtain  the  acceleration  components 
along  the  true  axes  of  the  accelerometer,  neglect  second  order  errors.  These  components 


are  then  substituted  in  the  model  equation,  Eq.  (2-1).  Using  the  criteria  of  Eqs.  (5-1) 
and  (5-2),  the  output  equations  are: 

a2  r 

A7r  =  H0  k  H1  acr  +  #  [  <K2i  *  Kt> <*  +  2<V  ‘  3K3i  »p  +  Sp  <*  '  2  V  +  3K3p  «l 

-(Kpo4V(0l-V*'/i’^  +  Kip]+(KM  +  K3p>  Jjf  <9*3> 

4s  f 

A8S  *  >0  *  h  %s  *  -f  1  <K21  -  Kt>  <>  +  2%»- 3K3i  Op  *  K2p  <*  '  2  0„> 

+  3K3p  »i  -  <Kpo  +  Koi>  <°i  -  “pi  V*  : k2  +  Kip  ]  -  +  K3p>  W 

Equations  (9-3)  and  (9-4)  show  that,  due  to  the  magnitude  of  the  allowable 
misalignment  angles,  the  value  of  KJp  can  be  determined  from  the  quadratic  terms 
only  if  the  magnitudes  of  the  products  of  the  accelerometer  coefficients  by  the  mis¬ 
alignment  angles  are  small  compared  to  the  magnitude  of  K.p.  It  may  be  found 
necessary  to  tighten  the  tolerance  on  the  misalignment  angles.  It  should  be  recalled 
at  this  point  that  the  misalignment  angles  Oj  (j  =  i,  p,  o)  are  unknown  and,  in  general, 
are  different  in  each  mounting  position. 

If  the  above  criteria  are  satisfied,  then-Eqs.  (9-3)  and  (9-4)  may  be  simplified 


A7r  H0  +  H1  acr  +  H2  acr  +  K3Ip  acr 
A8s  =  *0  “  X1  acs  +  l2  acs  “  K3ip  acs 


(9-6) 


where 


H2  -  f  [K2i  +  Kt+K2p  +  /ikz+  Ki 
I„  -  A  lK2i-Kt+K2p''/Tkz  +  Ki 


(9-7) 


K3(p=  <K31  +  K3p>/S'/T’ 
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■  ;v;':'- '' V>1^_  -  . -V 


The  coefficients  H^,  Ig,  Hj,  etc. ,  and  their  uncertainties  are  determined 
by  the  methods  outlined  in  Section  G.  The  value  of  K.p  is  obtained  from  Eqs.  (9-7) 


and  (9-8)  and  the  previously  determined  values  of  K_,  and  K  i.e. 
Kip  =  H2  +  l2  “  (K2i  +  K2p} 


(9-9) 


Even  though  the  parameters  of  the  output  equations  and  their  uncertainties  are 
such  that  Kj  cannot  be  determined  from  these  tests,  nevertheless,  the  centrifuge 
test  may  be  useful  for  determining  if  there  aro  hackling  ot-  ether  unexpected  effects. 
As  usual,  the  .residuals  should  be  plotted  to  see  if  they  are  randomly  or  systematically 
distributed. 


SECTION  X 


CROSS  PRODUCT  COEFFICIENT  K 

po 


The  best  mounting  positions  for  obtaining  the  cross  product  coefficient  are 
those  in  which  the  acceleration  components  along  PRA  and  ORA  are  reversed  in  sign. 
This  is  necessary  in  order  to  have  a  common  cubic  term  for  optimum  separation  of  the 
quadratic  and  cub.  terms.  The  compliance  term  in  the  quadratic  coefficient  of  the 
output  equations  are  reversed  in  the  two  positions  relative  to  the  cross  coupling  term. 

Positions  9  and  10  of  Fig.  4  satisfy  the  above  criteria;  The  acceleration 
components  along  the  reference  axes  are: 

Ai9  =  aw  “  1 "  'V  acr  1 


V  =  ‘  <au  -  V  =  y|-  [4  kz  *  acr  <l  +  k*  +  V  ] 
V •  (au  +  V  “  [4  %  'hr  b  V j 


<10-1) 


i_  “aw"  "kyacs) 


A_1ft  =  —  (a, +  a„)  =  Jk  I  K,  “  3.  (1 

plO  '  u  v7  cs  z  cs ' 


•vv  ] 


Aol0  =  Jz  (au"  V~  ~h  [  acs  kz  +  acs  (1  +  kx  +  ky 


(10-2) 


Substitute  Eqs.  (10-1)  and  (10-2)  in  Eqs.  (4-9)  to  obtain  the  acceleration 
components  along  the  true  axes  for  each  position.  These  components  are  then 
substituted  in  Eq.  (2-1)  in  accordance  with  the  criteria  of  Eqs.  (5-1)  and  (5-2). 
The  output  equations  are: 


Preceding  page  Wank 
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V'  J0  +  J1  acr  +  "F  [K2p  '  1  *2oi>  *  3K3p°o  *  K20  <l  '2  °i> 

1  “fcr 

+  3K3o  »p  *  'Klp  +  Koi>  Pp  - °o>  -  2ky  +  Kpo  |  +  <K3p  +  K30>  Jyj-  '1#-3> 
“cs  ( 

A10s=  L0  '  ll  “cs  *  -f  lK2p(1  +  2tV 


*  3K3p  %  +  K2o  <>  -  ^  -  3K3o4p 


-  <Klp  +  Kcl>  <“p  *  V  +  2  ky  +  K 


'poj 


cs 


‘■V'W  27j 


(10-4) 


Because  of  the  magnitude  of  the  allowable  misalignment  angles,  the  value  of 

in  the  quadratic  terms  of  Eqs.  (10-3)  and  (10-4)  can  be  determined  only  if  the  magnitudes 

of  the  products  of  the  accelerometer  coefficients  by  the  allowable  misalignment  angles 

are  small  compared  to  the  magnitude  of  K  .  As  previously  mentioned,  it  may  be 

po 

necessary'  to  tighten  the  tolerance  on  the  misalignment  angles. 


If  the  above  criteria  are  satisfied,  then  Eqs.  (10-3)  and  (10-4)  may  be 
simplified  to: 

A9r  =  J0  +  J1  acr  +  J2  acr  +  K3po  acr 
A10if  L0  “  L1  acs  +  L2  acs  “  K3po  acs 


(10-5) 

(10-6) 


where 


J2  '  f 


K2p  +  K2o  -  2  ky  +  Kpo 


L2  -  f  [K2p  +  K2o  +  2ky  +  Kpo 
K3po=  (K3P  '  K30> 1  2 


(10-7) 

(10-8) 


The  coefficients  J0,  LQ,  Jj,  etc. ,  and  their  uncertainties  are  determined  by 
the  methods  outlined  in  Section  6.  The  value  of  K  is  obtained  from  Eqs.  (10-7)  and 

pu 

(10-8)  and  the  previously  determined  values  of  Kgp  and  K2(),  i.e. 

Kpo“  J2  +  L2  -  <K2P  +  V  <W 


56 


A12»  *  No  -  N1  acs  *  2  [<K2l  '  Kt>  <>  -  2  V ''  3K3i  *  o  +  K20  <l  +  2  V  +  3K3o  »i 

*3 

-  <%  *  V  <°o  * °1>  -  &  \  +  Kolj  4  *  <K3o  +  K3!>  ^  (U-4) 

As  in  the  determination  of  K.p  and  ,  the  products  of  the  accelerometer 
coefficients  by  the  allowable  misalignment  angles  in  the  quadratic  terms  of  Eqs.  (11-3) 
and  (11-4)  must  be  small  compared  to  the  magnitude  KQi .  If  these  products  are  small, 
then  Eqs.  (11-3)  and  (11-4)  may  be  simplified  to: 


Allr-M0  +  Ml»cr  +  M2  4+K3ol4 

(11-5) 

A128=  N0  "  N1  acs  1  N2  acs  "  K3ol  aca 

(11-6) 

where 

M2  ’  i[K2l  +  kt  *  K2o  *  *  Koi]= 

(11-7) 

N2  '  \  [K2i  -  Kt  +  K2o  -  ^  +  Kol] 

(11-8) 

K3or  <So  +  K3i>/2^ 

The  coefficients  MQ,  NQ,  M^,  etc,,  and  their  uncertainties  are  determined 
by  the  method  erf  least  squares  as  outlined  in  Section  6,  The  valueof  KQi  is  obtained 
from  Eqs.  (11-7),  (11-8)  and  the  previously  determined  values  of  and  K2q,  i.e. 


Koi  *  M2  +  N2  -  <K2i  +  K2o> 
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